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>- ■ Abstract 

t^^ , Using a proper gauge condition the static spherically symmetric solutions 

^^ ' of Einstein-Maxwell equations with charged point source at the center are 

f— N . derived. It is shown that the solutions of the field equations are a three- 

■^ I parameter family depending on the Keplerian mass M, the charge Q and 

the bare mass Mq. The result can be interpreted as a correction to New- 
ton's gravitational potential and Coulomb's electric potential which are 
, both regular at the centre where the massive point is placed. A correc- 

g_) I tion to Gauss theorem is derived based on the nontrivial topology of the 

r^ • corresponding spacetime. 

> 
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-^ ! 1 Introduction 

Since the beginning of the 20th century just after the proposal of General Rel- 
ativity (GR) there were several attempts to describe the gravitational field of a 
physical point mass. It was a problem with interesting consequences in Maxwell 
electrodynamics and Newtonian gravity. The problem was that Maxwell electro- 
dynamics implies relativity principles while in case of Newton gravity they did 
not hold. After GR was invented the simplest solutions considering spherical 
symmetry and vacuum conditions were discovered. These are directly connected 
to the physical point source outside the source due to symmetry reasons but right 
after their discovery a new problem appeared. The vacuum solution of GR did 
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not fit into the new for tliat time matfiematical description of a physical point 
source. The latter was worked out by Paul Dirac who proposed a new class of 
unusual "functions" describing the point source and obtained as a limit of some 
real density distributions. They were called distributions (generalized functions) 
and have quite nontrivial mathematical properties PP, [2], [3j. These functions 
are strongly connected to the physical definition of a point source which is quite 
different from the mathematical one. Thus in order to understand the point 
source problem its the very definition that has to be considered. 

Coming from the time of Newton there is a significant difference between the 
meaning of geometrical point in mathematics and point mass or point charge in 
physics. First of all the abstract mathematical definition of a geometrical point 
is related to some dimensional properties, measure, etc. of this object when 
considered in the presence of other points, i.e. it is a zero dimensional object 
satisfying some extra assumptions when put into a set of other geometrical points. 
The corresponding physical definition is rather different instead. It rises when 
considering a finite object with well defined dimensions, for example it is a sphere 
of radius a, a cube with side length a or whatever else. This object is considered 
as a physical point only if an observer is placed in the picture at a distance L and 
the condition L ^ a holds. This definition strongly depends on whether such 
condition holds and leads to some interesting effects when point source solutions 
are considered. For example, in classical electrodynamics when static solutions 
are taken into account, the electric field ip of such a physical point charge is given 
by the simple relation ip ~ 1/L. If one tries to use this solution when the point 
source condition is not valid it is obvious that there is a significant divergence of 
the field. The question rising at this point is whether such solutions can be used 
when L ~ a or even when L -C a. The answer is well known since the invention 
of quantum electrodynamics (QED) - NO. It is because in such cases the physical 
structure of the "point" starts playing a significant role and one has to impose 
some quantum corrections to the electric potential. 

In case of GR however there were many attempts to derive a description of 
such a point source but none of them was satisfactory. In a recent paper |lj a 
way of finding such solution was described. It differs from previous attempts in 
the assumption that the well known gauge freedom of GR was used to regularize 
Einstein field equations. In many of the previous solutions this gauge freedom 
was neglected due to coordinate independence of all physical theories. However 
it is a well known mathematical fact that the use of singular transformations 
may lead to contradiction in the theory and thus has to be avoided. This is the 
reason for choosing only regular coordinate transformations in order to arrive at 
a correct theory. These regularity conditions in case of point mass lead to very 
interesting new effects. The main result is that the vacuum solution was cut 
to some value outside the coordinate dependant event horizon thus leading to a 
finite luminosity of the point mass in a very natural way. In other words this 
solution provided a regularization of the gravitational field of the physical point 



source and a corresponding correction to Newton's law typically acting only when 
approaching the gravitational (Schwarzschild) radius. 

Since the case of a point mass was regular it was interesting to try to solve the 
corresponding problem when the point is charged. Thus besides its gravitational 
field it produces also electric field in the static case. This article is devoted to 
derivation of such solution and the results are even more astonishing than the 
previous. First of all this method for the first time allows the description of real 
elementarr particles as well as astrophysical objects. Also besides the regular 
behaviour of the gravitational field it turns out that gravity regularizes in a 
natural way the electric field as well. There is also corrections to both Newton's 
and Coulomb's laws where the second is regularized at distances near the classical 
radius of the point source in case of elementary particle like solutions. Thus 
the interplay between gravity and electrodynamics opens new unexpected ways 
for regularization of classical divergencies and description of finite classical (not 
quantum) models of elementary particles. If this concept is right these classical 
corrections together with the proved to be finite quantum ones may lead to a new 
unified description of elementary particles. 

Having in mind the fundamental character of these considerations their de- 
tailed derivation and discussion will be delivered in the present article. 

2 Action Principle and Field Equations 

In order to obtain solutions of Einstein-Maxwell equations in case of a point 
source it is the spherical symmetry and the static character of spacetime that 
will be the starting point. It is well known that each symmetry of the action 
integral for a particular system reduces the degrees of freedom thus leading to a 
specific kind of dimensional reduction. 

In case of spherical symmetry in GR it is the group SO (3) that acts invari- 
antly as an isometry on the Riemann spacetime manifold 'M.^^'^'{gij{x)}. In jS] 
this invariance of the action integral was used to demonstrate that the latter 
can be reduced to two dimensional action represented as an integral with re- 
spect to the cosmic time t and radial coordinate r on the orbit space M^^'^^ = 
M[^^'^''{gij{x)}/S0{3). This reduction is a consequence only of the spherical sym- 
metry of the problem. 

Considering static solutions of Einstein-Maxwell equations in the rest frame 
of the source there is another symmetry of spacetime - it is invariant under one 
dimensional translations T(l) with respect to the cosmic time t. This additional 
symmetry reduces the problem to one dimensional so that it is described in terms 
of a single essential space-time coordinate - the radial variable r G [0, oo) as a 
proper coordinate on the one-dimensional manifold M*^^^ = Wl^^'^^T^l). 

The reduction of the classical action due to spherical symmetry is not the 
usual way to solve Einstein-Maxwell equations for a point particle but it will be 



the one used here. Such procedure was first introduced by Weyl and discussed 
by Levi-Civita jH] later. It yields the same solutions as the direct use of Einstein 
equations and is based on restriction of the classical action on the orbits of the 
5*0(3) group. The reason for choosing this type of description is that the action 
principle gives the simplest opportunity for solving the field equations correctly 
with their invariant distributions on the right hand side. These distributions arise 
from considering the particle's mass and charge density as proportional to the 
three dimensional invariant Dirac delta functions. 

2.1 Reduction of the Field's Action Integrals 

The most general form of the spherically symmetric line element is described by 
three independent functions in the metric tensor: 

ds^ = gu{r)df + g,r{r)dr^ - pirfdO^ - p{rf sin^ ed(j)\ (1) 

According to Birgkhoff's theorem [H], [7j the spacetime is static so that the metric 
components gtt > 0, g„ < and p > are unknown functions of the radial 
coordinate r only. 

Because of the gauge freedom in GR the gravitational field depends on only 
two of these components so that one of them can be fixed as a function of r. From 
geometrical point of view one has to introduce an additional mathematical struc- 
ture describing some specific embedding of the manifold M*^^-* into the spacetime 
iy[[(i,3)^ This can be done in different ways. If regular gauge transformations are 
used, the results do not depend on the choice of this embedding, i.e. the regular 
gauge freedom appears. However the use of singular gauge transformations can 
change the physical properties of the solutions, see |4J and the references therein. 
The most common gauge condition in solving Einstein equations is p{r) = r 
which was first used by Hilbert and provides the well known Reissner-Nordstrom 
form of the vacuum solution of Einstein-Maxwell equations. In the present ar- 
ticle this gauge freedom will be kept until a proper condition for gauge fixing is 
derived. The proper condition comes from regularity of the field equations, i.e. 
finite multiplier of the delta functions. 

2.1.1 Gravitational Field Lagrangian 

In order to reduce the action for the gravitational field without fixing the gauge 
proper second derivatives have to be extracted from the Ricci scalar. Then the 
Einstein-Hilbert part of the action, according to |lj, [H], gives: 

Agb^ = T^ / d^xJJ^lR = [ dtdrCcR, 
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where Cgr is the reduced lagrangian depending on gttiJ^), gwij') and p{r-) 
1 
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2.1.2 Electromagnetic Field Lagrangian 



Let's now consider the particle's rest frame so that the electromagnetic potential 
and the non vanishing components of the electromagnetic field tensor are 

A^{x) = (v9(r),0), J^rt = -^tr = v'{r). 

This gives a reduced action for the electromagnetic field 

Aem = —;]- f d''xJ\^\J^^,J^^''= [ dtdrCEM, 



H,r 



where the reduced lagrangian is given by 

^i.M = ^^=44L==[v.'(r)]l (3) 

2.2 Reduction of the Source's Action Integrals 

In order to obtain correct solution, the corresponding contributions of the point 
particle's self energy and the interaction with its own electromagnetic field have 
to be added to the total action. Later on it will be shown that the second, the 
so called electromagnetic self interaction, is well defined because of the regular 
behaviour of the electromagnetic field at the center of symmetry r = 0. In other 
words our solution of this problem shows that there are no divergencies of the 
electromagnetic field at classical level. This is a consequence of the nontrivial 
geometry of spacetime. This may turn to play a role of a universal regularization 
factor because an infinite density of energy of any kind according to GR will 
cause similar deformations of spacetime. It turns out that this phenomena comes 
with another one which is disappearing of the horizons, i.e. no divergencies in 
the gravitational field either. These are the most important new effects under 
consideration in the present article. 

The first of the above mentioned action integrals is well known and it can be 
written in the following form: 

Amo = - e{r)\[\^\<fx = -Mq I ds= dtdrCu^- 



t.r 



Here L*^^) is the one-dimensional world-line of the particle and the second relation 
comes from the use of the three dimensional invariant Dirac delta function for 



the particle's mass density: 

e r = Mo6Jr) = -7^= = J— \ no - (4) 

The mass Mq in this formulation is the "bare" mass of the source corresponding 
to the energy of the static configuration. In other words it determines the energy 
needed to form the point source. Having this in mind the reduced lagrangian for 
the self energy is 



Cm, = -Mo^/^^Sir). (5) 

The same consideration can be used to reduce the action describing the self 
interaction of the source 

Aq = - [ d^x\/i^\A^ix)fix) = -Q f A^u^ds= [ dtdrCq. 

This action integral and the action for the pure electromagnetic field are both 
manifestly covariant under diffeomorfisms of spacetime and U{1) gauge group in 
electrodynamics. This is an important feature because it provides invariance of 
the solutions under such transformations with one restriction - to consider only 
regular gauge transformations. 

The charge current and charge density are given by 

f = gir)u^, g(r) = Q6,ir) = -^^ = --2!^ , 

where the velocity u^ is a unit length four-vector. If the particle's rest frame 
is chosen then u'^ = {\^g^, 0) so that the reduced lagrangian describing the self 
interaction is 

Cq = -Qv{r)5{r). (6) 

2.3 Field Equations and Quasinormal Coordinates 

To derive equations for the components of the metric tensor and the electromag- 
netic field of a point source the well known action principle will be used. This 
means that one has to sum equations 0, ©, © and (jH)) and find the variation 
of the reduced total action integral: Atot = Jm^^-^'> ^^ drCtot where 

Ctot '■= Cgr + Cem + Cmo + Cq. (7) 

Before doing this it is useful to transform the field variables gtt{r), grrir) and 
p{r) to three new functions of the radial coordinate - V9i(r), ^2{f^) and <^{r): 
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v^tt = e''^ v^^^ = e^^^-^1-^ p=-pe 

Lfi = In(v^), (/?2 = In ^-^ , ¥? = In _ 

\ P J VP v^^fc/ 



(8) 



which partly diagonahze the total lagrangian (|7j) and will be referred to as quasi- 
normal coordinates. The symbol p denotes a normalization constant to be defined 
later from convenience reasons. The total lagrangian in terms of these new func- 
tions has the following structure: 



r -1 

'~'total „ 



p^e^ ( - if'^ + if'^ + e-2^V") + e'^e''^' - (MqC^^ + Qip^ 6{r 



This lagrangian does not depend on derivatives of (^(r) so that the factors e^'^*-'"-' 
can be treated as Lagrange multipliers. This result comes from the above men- 
tioned gauge freedom that was left. It can be fixed imposing a constraint in the 
lagrangian. According to the theory of constraints the field equation obtained 
from variation with respect to cp define a constraint on the solutions. 

The full set of field equations are derived from variation of the total action 
with respect to all fields - ipi, ip2, (p, and p: 



dr \ dr J p^ 

dr \ dr J p^ 



fe\- \p'l + \p'^ + \e-'-^p'') - ^e^-e- ^ 0, 

where the sign = was used to denote the constraint. 

Finally changing the electromagnetic field p{r) to a new variable v'3(r) : 
p'^{r) = (y9'(r)e~'^^^''^ the full system of differential equations for the new fields in 
arbitrary gauge are 

Arp^=p',' + ^e^^~n{r), 
P 
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A,<^2 = ^e^'^^-^^ (9) 

P 

A^V^g = p[p'^ - —e'^'~'^5{r), 

P 
1 / \ 1 

2' 



;pV(-^f + ^? + ^f)-V^e-^ 



In these equations the differential operator A,. := e~'^j^ (^'^^) ^^^ introduced. 
This operator is connected to the radial part of the three dimensional Laplace- 
Beltrami operator: A,. = —g„Ar, where A^ comes from the definition of the 
laplacian in curved 3D spacelike slices at a fixed time t: 



A = -l/^\^\d, [yWW'd,) = -g'^ {df,-f,d,) = Ar + l/p'Ae^. The four 

functions T^ appearing in this expression are not components of a vector field 
and their definition according to [S], is 

One can fix their local values in an arbitrary way thus fixing the gauge of the 
problem, i.e. spacetime coordinates. In case of static spherically symmetric met- 
ric these functions are {f t, f r, Tg, f ^j = {0, (f', ta.n9, 0} which again shows that 
the variable ip{r) is the one fixing the radial gauge freedom. The global conse- 
quences of such gauge fixing need a very careful study because of the distributions 
appearing in the field equations and the possible divergencies they can provide. 

3 Gauge Fixing and First Integrals of the Field 
Equations 

In order to have regular field equations one has to find a proper gauge condi- 
tion. The proper choice of (^(r) must insure a finite multiplicative factor e'^'^^^^~'^^'^^ 
of the Dirac delta function when r goes to zero: 

lim |(y9i(r) — (^(r)| < oo. (10) 

r— >0 

This provides regular field equations at the place of the source r = ^. The 
easiest way to find a proper gauge satisfying this is to solve the system (jH} outside 
the source and take the limit r —>■ with some assumptions for the asymptotic 
behaviour of <^(r). 

3.1 Solutions Outside the Source 

The first step in deriving solutions outside the source is to change the radial 
coordinate r to a new variable a : da = e~'^dr. This transforms A^ to a new 
differential operator acting on functions of the new variable a: e"^"^^ so that 

"'^It is important to note that the center of spherical symmetry r = is mathematically 
excluded because of a coordinate singularity. In many problems it can be added back to the 
manifold by a trivial mapping but this problem requires a more sophisticated procedure due to 
the 3D invariant Dirac delta function placed there. 



the field equations Q outside the source in terms of a read 




(11) 



0. 

Next it is useful to define variables which can be interpreted as canonical 
energy for the corresponding field such that 



A proper substitution of these in (jllll combines e\, 62 and £3 into two first in- 
tegrals: £2 = ^2 = const and ei + 63 = —€2- As usual these two first integrals 
are used to divide (fTT|) to two separate equations which in the current case are 
exactly the same while the equation for ip^ still depends on ipi: 

d'^ ,0 2e2 ,„ f dip-i\ / duDi\ 269 



da'^ p^ ' \ da J \ da J p^ 

These are the final field equations for the charged point particle outside the 
source. Depending on the sign of €2 this system has three different types of 
solutions. Later on it will become clear that this sign corresponds to charges 
higher, equal or lower than the mass in terms of the usual parameters used for 
description of the solution. These three solutions are: 
I. In case 62 > 0: 

^lia) = Ci — In I sinh( — - — a + C2 
\ P 

V^2(o") = C3 — \n i sinh( — - — a + C4 
V P 



II. In case €2 = 0: 



tanh(^a + C2) 



ipiia) = -\n{Cia + C2), 
ip2ia) = - In (C30- + C4) , 

<^(^) = T ^^ ' +C,. 
Gicr + G2 



III. In case €2 < 0: 



(flier) = Ci — In ( sin( — a + C2' 

P 



(^2(0-) = C3 - In ( sin( — ^cr + d] 



^{a) = T- 



oC\ 



tani 



v^=^ 



a + C2 



+ C5 



3.2 Basic Regular Gauge for the Point Particle 

After solving © outside the source one has to find ip{r) satisfying (fTn|) . A large 
class of gauge conditions is described by lagrange multipliers e^'^'-^-' with the 



following behaviour when r — > 0: e 



-ip(r) 



kr"". Considering this gauge class 



which is general enough there are two possibilities for the power n of r in the 
asymptotic behaviour 

I) If n 7^ —1 is chosen it implies a = -^^r"^^ so that \(pi{r) — <^(r)| is 



Ci + In A; + In 



sinh(^^r"+i + C2) 



In /c + In 



n+l 

Ci + In A; + In 



fc^^n+l + C2 



sm 



■fcv^ 



(n+l)p 



^r"+i + C2) 



for each of £2 > 0, £2 = 0, £2 < respectively. This means that the only choice 
for n that satisfies the gauge condition p()|l independently from the integrating 
constants is n = 0. 

II) The case n = —1 corresponds to a = fclnr which always gives divergent 
value of \(pi{r) — <^(r)| when r goes to zero. 

Finally the only choice that keeps the fields regular at r = is e"*^^*"^ ~ const 
for r —>■ 0. This defines a whole class of possible gauges which will be referred as 
regular gauges for the point particle source. The easiest way to choose one from 
those regular gauges is to impose the condition e""^^^^ = 1 or simply 



(f{r) = 0. 



(13) 



This choice is the so called Basic Regular Gauge (BRG) |1]. It will be used from 
now on for solving the Einstein-Maxwell point particle problem. 
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3.3 Field Equations in BRG and First Integrals 

The next step after fixing the gauge freedom according to (|T!?|l is to solve the field 
equations (jH)) in BRG. The latter read 

^2 = -,e'^' (14) 

P 

^3 = ^W, - ^e-^5(r), 



^p^(-^f + ^? + ^?)-V^ = o. 



In order to divide these equations into independent parts conserved quantities 
have to be extracted. These are first integrals for this system and they can be 
easily derived from combining (|T^ where the choice of BRG gives a = r. In 
terms of £1,2,3 (fT^ becomes 

^2 - ^'[ - K - ^'D = -I (^1 +e2 + £3) - ^e-^5(r), 

\ / pZ pZ 

4 = 0, (15) 

(£1 + £3)' = -4 (^o<^'i + Q<^3) e'^'Sir), 
P 

ei+e2 + es = 0. 

Because of the Dirac delta functions in the latter it is useful to introduce a new 
quantity A(r) := 6{r) — 1 which has the following properties 

1) A(")(r) = 5("-^)(r), supp(A(")(r)) =0, n>l 

2) (A(r))" = (-l)"-iA(r), n > 1. (16) 

This is a function which is zero for all values of r > and —1 for r < 0. 

With the use of (jl6|) and the definition of quasinormal coordinates (jH)) the 
equations (fT3j) transform to 

VaMgttir) p j p ' 

£tot + 4 (^o¥^'i(0) + Qv'M) e'^^^°)A(r) = E^ot, 

^tot '■— £1 + £^2 + £3 — U, 
^2 = £2, 
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where £2 and i?tot are constants. A combination of the second and third of these 
equations imphes 

so that the choice of a domain for r such that r G [0, 00) determines only one 
possible value of Etot = 0. With this value and the fact that A(r) = 0, Vr G 
[0, 00)^ there are three first integrals for the system of field equations 

— C = const, £2 = ^2 = const, £1 + 62 + ^3 = 0. (17) 



^/gJjJgJj) 



Finally using the second and third of these and the definition ()12|1 the system 
transforms to a very "easy to integrate" one: 



<*%-.. _^e-.-^«W = 0, 



dr"^ p2 


P^ 


d' g-,. 262 , 

dr^ p^ 


= 0, 


J2 12 2e2 
¥^3 =^l -2 • 





;i8) 



4 Solutions for the Gravitational Field 

The solutions of the field equations corresponding to point particle are these 
solutions of (fTHj) that satisfy the proper asymptotic conditions, i.e. the physics 
of the problem will determine the integration constants. It was shown that the 
general solution of the field equations (fH^jl outside the source depends on the sign 
of e2 so that there are three cases that should be considered: 

I) Positive values, 62 > 0. 

Later on it will be shown that these values of e^ correspond to charges smaller 
than the mass of the point source. The solutions of (P^ for the gravitational 
field in this case are 

g-^iM _ Q^ gj^i^ m (ri _ r)) + ^ sinh (fir) A(r), 

p^\l 

^-^2{r) = 6-2 sinh (fi(r2-r)), 

where fi := ^/%- and Ci, C2, Vi and r^ are constants coming from integration of 

the equations. 

II) Zero value, €2 = 0. 

^The function is zero but its derivative is the Dirac delta function 

12 



This corresponds to equal mass and charge and the solutions are 

^-Mr) = Ci (n -r) + ^rA(r), 

g-c^2(r) _ ^^ (^^ _ ^^ _ 

III) Negative values, 62 < 0. 
This is the most physically significant case because it corresponds to charges 
greater than the mass of the point particle. This is just the observed relation for 
real world particles. Here the solutions for the gravitational field are 

g-c^iM _ (j^ gj^ r^ ^^^ _r)] + ^ sin (fir) A(r), 



-^2(r) = C2sm[n{r2-r)], 



where fi := J-^- 

The integration in each of these cases is not straightforward because of the 
delta function. The way to solve (|18|) is to fourier transform the equation with 
a delta function and find a proper solution for each case. The final solution is 
derived by adding this result to the solution of the homogenous equation. 

4.1 Asymptotic Behaviour 

As was mentioned above the integration constants should be obtained from some 
physical conditions on the general solution. The first two such conditions are 
asymptotic fiatness of the metric and coincidence with the vacuum solution in 
terms of t, p, 9, outside the source, i.e. coincidence with the solution written 
in the usual Reissner-Nordstrom form. These are quite common while the third 
condition is less trivial and was recently discussed in |1]. 

The asymptotic behaviour of the solutions for the gravitational field is defined 
via the limit p{r) ^ 00 so that it should be obtained for each case individually. 

Let's first consider the case €2 > 0. Having in mind the expression for p{r) in 
terms of ipi and ip2 (|H1) one has 

Ci sinh [n (ri - r)] + ^ sinh (fir) A(r) 
C2 sinh [fi (r2 — r)] 



and it is obvious that p{r) -^ 00 when r — ^ r2 =: r, 
In the second case €2 = one arrives at 



oo- 



and again p{r) -^ 00 when r =— > r2 =: Too- 
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Finally when €2 < 

Ci sin [n (ri - r)] + ^ sin (fir) A(r) 
C2 sin [fi (r2 - r)] 

so that p(r) -^ 00 when r ^ r2 — ■^. In order to have divergence- free solutions 



one has to choose the value of k such that ^ — ^ ^ (^ — [^]) ~- '"oo' where 
means the integer part of -. This is so because using this choice r belongs to 
the interval (0,roo] and there are no divergent points in it, i.e. Q{r2 — r) G [0, vr) 
when r G (0, Too]- 

These results are quite significant because they mean that the physical infinity 
for the problem is at a finite value of the radial coordinate. In other words 
the regular gauge condition determines a finite interval (0, Too] for the radial 
coordinate in which the regular point-like solutions live. Later on it will become 
clear that this restriction for the domain leads to a finite value of p(r) when the 
center of symmetry r = is added to the manifold. This means nothing else 
but finite value of the luminosity of the point source. The reason to refer to p{r) 
as a luminosity variable is that it determines the luminosity via the expression 
L = 4^7-732 5 where I is the intensity of the source. 

At this point it is crucial to note that the finite domain for r is where arbi- 
trary regular coordinate transformations (gauges) are allowed. This is due to the 
regularity conditions that each gauge theory is forced to satisfy, i.e. regularity 
of the gauge transformation in the domain where the physical fields live. Thus 
up to now only the domain of the gravitational gauge invariant fields, the metric 
components, was defined. This domain is easy to be mapped on a new one, where 
r varies from to 00. This can be done by using a proper reg"u/ar transformation. 

After defining r^o as the end of the regular physical domain for r one has to 
impose flatness of the solution in the limit p{r) ^ cxd or equivalently r = r^o. 
This means that at r = Too the metric components have to coincide with the 
flat ones: 5'tt(^oo) = 1 and Qppir^) = — 1, where Qppp'^ = g„- Now using these 
conditions and (fTTj) one obtains 

- , "'"-) 

\/-9rr{r)gtt{r) 

Finally expressing gtx,{r) and the above first integral in terms of ^pi{r) and ^P2{'r') 
the asymptotic flatness gives the following two equations 



g-Vl{r-oo) ^l^ p 



(^-Vlij-W ^-V2{r) _ (^-V2{r)\' ^-Vi{r) 



Solving these determines the constants Ci and C2 for each of the three cases and 
the solutions of (|T8|) satisfying the above mentioned asymptotic conditions are: 
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1) e>0 



2) e = 



,„ t,.\ sinh [i7 (ri — r)l Mn , ,^ , ^ , . 



ri - Too p^ 

P 

3) e<0 

,„ / N sin [f] (ri — r)l Mn ,^ ■, ^ , . 

^ = • rm % + ^ ^^^ ^^ ^ ^ ' 

sm[i2(ri -roojj p^ii 

4.2 Coincidence with the Vacuum Solution Outside the 
Source 

The next step after determining Ci and C2 from asymptotic flatness is to impose 
that outside the source the solution for g^t is the same in terms of (t, p, 6*, 0) as the 
well known form of the vacuum Reissner-Nordstrom solution. One should note 
here that the coordinate transformation from the latter to the current solution is 
defined only in the domain r G (0, Tqo] which maps to p G (po, 00) with a non-zero 
value of Po = p(0). This means that the current solution is regular for r = and 
later on will be shown that it has no singularities in the region p & (po, 00). 

1) Consider the case €2 > 0, so that ^2 is defined as i7 = \/^- The solutions 

for gtt{r) and p{r) outside the source are 

sinh^ \Q{ri - r^o)] , , _ sinh \VL(ri - r)l Vtp 

9tt{r) = . , 2 r^, ^, P{r) = P- 



sinh^ [fi(^i — t)] ' sinh [r2(ri — Too)] sinh [(]( 

After simple algebraic calculations the following relation can be derived 

^ l]p^coth[^(ri-roo)] ^ (l]p^sinh-^[l](ri-roo)])' 

gtt{r) = l-2 '-— '- + y— ^. (19) 

p{r) p[rY 

Comparing this result to the well known expression g^^ = 1 — ^^ + % one ends 
up with the following system for the unknown parameters fi, ri, ri — Too and p 

fip2 coth [Vt{ri - Too)] = M, 



sinh \VL{ri — r 
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\Q\- 



Here M is the Keplerian mass fU] of the source, which in general is different 
from the "bare" mass Mq, and Q is the electric charge. The solutions of these 
two equations for Q and ri — Too are 



il = — a/1 - 772, coth [fi(ri - Too 
where 77 is a dimensionless parameter such that 



1 



q_^2 



V 



M 

M 



(20) 



Because of 62 > this solution restricts t] to values smaller than 1 which is simply 
the case of Q < M. 

With these parameters set the solutions for the gravitational field correspond- 
ing to ?7 < 1 are 

cosh [Q (roo — r)] H -j^=^= sinh [Q (too — r)] + 



-1^1 (r) 



+ 



Mn 



: sinh (fir) A(r), 



-1/32 (r) 



M^l-r/2 

sinh [fi (roo -'^)], 



A/ 



P 

2) In the case €2 = there is an expression for gtt{p) similar to (|T9|) is 



^tt(r) = l-2- 



+ 



(ri-rca)^ 



p{r) p{ry 



(21) 



which corresponds to M = Q so that ri— re 
of the field equations are 



M' 



According to this the solutions 



-ipi(r) 



-•fi2{r) 



1 + 



M 

7 



roo — T H rA r 



Tr^-r 



3) The case €2 < is very similar to 1). The difference is in Vt 
the equations for the integrating constants are 

fip^cot [fi(ri 

fip2 



ro 



M, 



sin [fi(ri - roo)] 
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— ^ and 



The solutions of these two equations are 



n = %V¥^, cot [n{n - Too)] = ^i=, (22) 



where the dimensionless parameter rj (J20|) is greater than one corresponding to 
e2<0. 

The final expression for the solutions is 

Q-viir) ^ cos [n (roo - r)] H sin [Q (rpp - r)] + 

+ sin (Or) A(r) 

_L 

g-v'zW _ ^^^sin[0(roo-r)], 



M 



^ := -^ 



4.3 Gravitational Mass Defect of a Point Source 

It is a well known fact that there is a difference between the mass in the Einstein 
equations and the so called "bare" mass which is the invariant quantity defined 
as an integral of the mass density. The inequality of these two masses defines a 
new free parameter in the solution which will be considered in this subsection. 

First of all in case of spherical symmetry one has the following expression for 
the bare mass 

Mo = 47r /" v/-fc(r)p2(r) e{r)dr, (23) 



where e(r) is the invariant density distribution defined in (@]). This is the quantity 
describing the energy of the particles configuration as discussed above. 

Besides the bare mass one can consider the mass in Einstein equations which 
is the so called Keplerian mass of the source and is the one that an observer at 
infinity will measure witnessing the motion of a test particle. In case of a point 
source the Keplerian mass is 



M = 47r / p\r)p\r)eir)dr = Mo / ^/^^S{r)dr = Mq^/^Jo). (24) 



The dimensionless quantity g := ^/gui^) takes its values in the interval [0, 1] 
and will be referred to as "gravitational mass defect ratio". Later on it will be 
shown that the values g = 0,1 are not physical, i.e. they give divergencies in the 
metric components. 
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Now using the definition of q and performing a rational change of the radial 
coordinate r -^ (^ which maps the finite domain r G (0, Tqo] to (^ G (0, oo): 

c - r^. ^ - -^ (25) 

gives the following solutions depending on sign(r7): 
1) ?7< 1 



^-'""-i7!=f-K^l' '''' 



T Ti A / 1 7^^ 

A$^ := arcsinh-^^ arcsinh-^^ 



2)r/ = l 



^^ 



-^i(c) . , 1 - ^ / 1 , 1 C 






3) r7> 1 



/A$+x sinfllf) sinfA<l>+^) 
e-MO = cos ( ^ ) + -ASlU + _L^1V A(C), 



eVT^T VC 

A(P := arcsm-5^ arcsm , 

1 

where ^ := ^ in all cases. 

The consideration for r^ > 1 is the most interesting one because it corresponds 
to charges higher than the mass of the point particle. This is the most common 
case in real world particle physics. In this solution there is a limitation for the 
charge of the point particle given by the last inequality. It comes from the relation 
Q'fiio) = \j Q between firoo and the dimensionless parameters: 

sin I arcsin — h r2roo 1 = — . (29) 



This equation has real solutions for Qvoo only if the right hand side is smaller 
than 1 so that a/^^ — 1 < f]Q- It is interesting to note here that this can be 
considered as a limit for the charge of the point particle: 

101 < " 



'I 



n2 



This limiting value depends on the gravitational mass defect g which is very small 
for real particles because of their weak gravitational field so that their charges 
can be big enough. 

The luminosity variable in all these cases does not depend on p and it is useful 
to define a dimensionless quantity ^(C): 

o(C) e-'^i('^) 



which is just p{() measured in units of M. This function outside the source is 
given by: 

. v<i, z{o = 1 + v^r^^coth (^), 

. ^>l,^(C) = l + v/;^^^Icot(f}f). 
The value of z{0) is the same for each of these three cases: 

2 

where in case 2) z{0) is given by the simple substitution r] = 1 in the last expres- 
sion. 

In case of particles, with r/ > 1 a limiting value of the gravitational mass 
defect ratio exists: 



which can be derived in the same way as the charge limit discussed above. The 
physical relevance of this limiting value becomes clear after its substitution in 
fl30|) . It gives a limit for the initial value of the luminosity variable which is 

pT = %■ (31) 

This is exactly the classical radius of the particle derived from elementary energy 
conservation assumptions - pci = Q"^ /M. Within the current framework this 
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radius corresponds to the minimal value po ^ Po™ of the luminosity variable 
due to the gravitational field produced by the point source. In the next section 
it will be shown that this provides a cut for the corresponding electromagnetic 
potential and thus regular behaviour at the point r = 0. This leads to a cut of 
the electromagnetic field if the motion of a charged test particle in the presence 
of point-like source is considered. 

Another interesting consequence when ?; < 1, is that the cut in p(r) is al- 
ways outside the event horizon appearing in the vacuum solution without any 
restrictions for the gravitational mass defect ratio g: 

z{0) > 1 + Vl-r/2. 

This corresponds to disappearing of the coordinate singularities which arise in the 
physical domain in the classical Reissner-Nordstrom form of the vacuum solution. 
For our solution the horizon and all phenomena related with it are in the non- 
physical domain of the variables. It is important at this point to note that the 
solution derived here and the Reissner-Nordstrom one are the same outside the 
source in terms of t, p, 9 and cf) but there is a significant difference between them 
because of the above mentioned cut on p{r) . In case of the classical solution there 
is a S{p)/p pathological divergence at the point source placed in p = while the 
current solution has only 6{r) type singularity. This is so because of the choice 
of regular gauge when describing the point source placed at r = 0. 

Finally one has to consider the limits on the gravitational mass defect. They 
can be derived from regularity restriction for the metric components. In case 
oi g = the multiplier of the delta function e"^^^' in (|T8|l is infinite while in 
case of ^ = 1 the cut p(0) becomes infinite. Both of these cases have to be 
excluded from the domain in which g varies because of these two singularities 
and then the restriction for g becomes stronger: g G (0,1). It is the case of 
strong gravitational fields that take place when g approaches 1 and the asymptotic 
behaviour ^ — ;> 1 ^ po = oo has a simple physical interpretation. It is connected 
to the luminosity of the point source in such case, which tends to zero due to 
L < Lo ~ Pq'^ and thus stronger gravitational field corresponds to darker object. 
This is exactly the same physical situation as in purely gravitational point source 
|1] where po = 2M/(1 — g"^) and the difference between the two values of po is 

1- /Einstein Einstein-Maxwell '\ ^ 

,™ ^^0 ^0 ^ ~ 2M' 

where Einstein stays for pure gravitational field and Einsten-Maxwel for the 
current solution. 

5 Solutions for the Electromagnetic Field 

The gravitational field of the point source is given in (j2ni), f)27p and ()28p with the 
use of (jHl). Now one has to consider the electromagnetic field which is a solution 

20 



of the third equation in ()14|) with ^'^{r) = Lp'{r)e 'f^^'^\ It is easy to check that 
this equation is equivalent to 



y,'3(r) = y?'3(0)e^^('')-^^(°) - ■§e^^('-)A 



P' 



r). 



Going back to the electric field '^{r) this is transformed to 

^'M=(«AW-p-^)(i)' (32) 

and straightforward integration gives the solution for ^p{r) depending on the initial 
values ipifi) and (/''(O) 

^W = ^(0) + Q{-^)-ik?l ^(^) - p'^-^ {w)-m)- (33) 

To find the initial values one has to compare ()33|) to the well known static 
electromagnetic field outside the source, i.e. for r > 0: 

V?(0) + 



p(0)f?^ / p{i^)q'^ p(^) 

This is equivalent to two equations for the constants ^^(O) and v^'(O) and after 
solving them and performing the rational transformation (j^H|) the final expression 
for the electromagnetic field becomes 

The most significant point about this solution is its finite value for C = 



^(0) = sign(Q) ^-V'-"''^-^' . (35) 

V 

This shows how an appropriate regularization of the gravitational field of a point 
particle suppresses the divergencies of the electric field as well. This purely GR 
phenomena is a regularization at classical level which is impossible to derive in 
terms of Newtonian theory of gravitation. The importance of this result comes 
from the independence of this regularization scheme on Maxwell theory. This 
gives a chance for regular solutions of the point particle problem in generalized 
theories coupling these two like 5D Kaluza Klein, extra dimensions, etc. 

Another interesting feature of this cut for the electromagnetic field comes 
from the restriction on the luminosity variable when the charge is greater than 
the mass of the source, i.e. when real particles are considered. In this case it was 
shown in the previous section that there is a limiting value for the mass defect 
parameter g^^ and thus p(0) > pd which corresponds to ip{0) < Q/pd- Thus the 
electromagnetic field of a test particle with the same charge as the point source 
is restricted to values smaller than the Keplerian mass of the source M. 

21 



6 The Final Solution 

In order to arrive at the final form of the point particle solution one has to perform 
a coordinate transformation of the dimensionless quantity ( defined in fl25p to a 
new radial coordinate. The easiest way to do this is by introducing a constant f 
such that 

r 

This is a rational transformation which is well defined in the domain ( G (0, cxd) 
and will not lead to any poles in the solutions even when the point ^ = is added. 
The scaling parameter f should be such that gtt has the standard asymptotic 
behaviour when r -^ oo: gftt ~ 1 — — . This condition gives different expressions 
for f in each case rj < 1, t] = 1, t] > 1: 

/l 2 / 2 1~ 

r = ^ ~^ M, 77 < 1: f = -^—M, ?? = 1: f = ^^ ~ M, ri>l. 

Using this result and dr = ttty^' coming from (|^H|) the point particle solution of 
Einstein-Maxwell equations can be written in the following way 

ds^ = e^^'^W (dt^ - -j^) - p{r)^ [dO^ + sin ^2^02^ , (36) 

Q 



(p{r) 



p{r) 



where the notation {pa{r; M,Mo,Q) = (pi{r) was used. Each case for r] leads to 
different expressions for each of the functions (p{r) and p{r) and there are three 
different solutions determined by the following functions: 
1) In case M > Q 



(fa{r) = In 



r+fi r + ri 



^si„h (4p! "f^' , 



(37) 



sinh 



y/M^ 



pW) = M + JM^ - g2 coth I ^^^ — ^)+Mo ^,- '^'^ / A(r). 

V ^ + ^^ ) sinhf^^ 

where the scaling parameter fi depends on M, Mq and Q: 

v/M2 - Q2 



ri : = 



arcsinh ( — — ^ ^ ) — arcsinh ' — 
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2) In case Q = M 



p{r) 



In 



p{r) 



M + f2 + r + Mo— A(r) 
M2 



(3^ 



r2 



Mo-M 



3) In case M <Q 

( 

(fc^ir) = In 






\ 



r + r3 



V 



sm 1 -f^' \ P^'^ 



(39) 



/ 



p(r) 



M+v/Q2T:M2 cot 



where a different scaling parameter fs is used: 



sm 



Mn 






sm 






Mr), 



n ■-- 



VQ' - M2 



arcsm 






Mo 



— arcsm 



Q 



In this case there is also a restriction for the charge of the point particle: 

M 



Q< 



/mV 

\Mo) 



providing that the minimal value of p{r) is the particle's classical radius. 
The function Na{r) in each case is given by 



Nair) 



M {r + r) 
r p{r) 



(40) 



in terms of p{r) and f = 'ri,2,3- 

All these functions are defined in the interval r G (0, cxo) and are singularity 
free everywhere in this domain. A significant property of the considered solution 
is its regularity even when r ^ which is very uncommon due to the well known 
divergences emerging in the classical theory. 

The functions v^a(r) and ^p{r) with the use of p{r) in each case correspond 
to corrected Newton and Coulomb potentials of the point source respectively. 
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In order to investigate the obtained corrections it is useful to rewrite the scahng 
parameter for each case in terms of po = p(0) = Mz{0) (see (j3Uj) for the definition 
of z(0)): 



n = — , ^ , r2 = Po- M, rs- 



arctanh-^^ irr- arctan — -rj— 

Now using these one can consider the following asymptotic cases: 

• Vanishing electric charge Q ^ and thus Q <^ M. This case leads to the 
pure gravitational point source solution derived in |3]: 

2M M 



Po -^ -. ^ =^ ri 



1-^2 InMo/M' 

and the corresponding corrected Newton and Coulomb potentials are 

<^"(0 , M ^ <^(0 ~ 1 - e ^^^wi^ ^ (41) 

^ + InMo/M V / 

• Mass of the point source M <^ Q where the condition r/"^ + ^^ > 1 is 
satisfied. This case leads to the following asymptotic behaviour of the 
scaling parameter -r^: 

where the limit is considered over the the curve 77"^ + ^^ = 1. The corre- 
sponding corrected Newton and Coulomb potentials are 

(fair) ^ - In cos ^ J , (p{r) ^ tan ^ (42) 

7 Hamiltonian and Energy of the Fields in BRG 

After solving the field equations using the lagrangian formalism one can step 
back to the beginning and derive a hamiltonian density corresponding to ((Zj). 
In the scope of the present article this will be used to find the total energy 
of the charged point source. This procedure was used in |lj in case of neutral 
point source and will be repeated here. If BRG is taken into account the total 
lagrangian ((7j) depends on three "dynamical" fields ipi, (p2, ^3 and one can obtain 
their conjugate momenta from ttj := ^^f : 

TTi = -pV'i, ^2 = PV2, ^3 = p'e-2^V'- (43) 
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Using these relations the total hamiltonian is simply 

'Htot = ^ vTiV^i - £tot(vr, if) = etot + (Mdc'^i + Qip) S{r), 

i 

where fll2|l was used. The total energy is an integral of this hamiltonian density 
over the whole domain of the radial coordinate r: 

^tot := / nUr)dr = M + M 1 - ^1 - (1 - g^W 
Jo '- 

The second part of this expression comes from the fact that etot = ^i + ^2 + ^3 
satisfies the third relation in (fTTj) . 

Next it is useful to separate the energies corresponding to electromagnetic 
and gravitational field. The easiest way to do this is using the energy density of 
the electromagnetic field derived from its energy-momentum tensor: 

cem = Temo = ^r;^- (^^) 

The electromagnetic part of the total energy is then determined as 

^EM = / ^TiphEudp = ^ = M 1 - ^1- (1-^2)^2 _ (45) 

Jpo PO ^ / 

The energy of the gravitational field can be obtained from the relation jl] 
Etot = Mo + -Egr + Eem so that 

Egr :=M-Mo = -(l-f?)Mo<0. 

This result was recently obtained for the pure gravitational field of a point source 
without electric charge in ^j and it turns out that the charge does not change it 
at all.. 



8 Integral Theorems and Invariant Quantities 

After fixing the gauge condition ()13|) the two field equations containing delta 
functions are 

P 
Ar^, = ^[^',-^e^^5ir). (46) 

These can be used to derive integral invariants for the fields just like the Gauss 
law in case of electromagnetic field only. Actually one of these integral invariants 

25 



will be just the Gauss law with a correction due to General Relativity. In order 
to derive these invariants one has to transform (j46j) to equations for ipi and if 
corresponding to the gravitational and electromagnetic field respectively. After 
having this done the equations read 



Aripi 



Arf 



^ 



/2 



Mo y/gTt , . 



-9rrgtt P" -gvr 



-9rry/9tt 



P^ -9tt 



Sir). 



A combination of ()32|) and (fTB|) as well as (jH)) in BRG give if' = —Q9{r 

2 , 

and _2 )— ^ = 1 and using these the field equations are transformed to the fol- 
lowing 



Arip = 

Due to r G [0, oo) this reads 



1 



gtt + MoS{r) 



VW\ 



Q^ejr) 

p' 



[2{VKt)'0{r)-^,6{r)]Q. 



4:'K\ ng^^^Aripi = Air 



in 



Q' 



i-Kp" 



\[\g^\+Mo6{r) 



A'KJ\gi^)\A.,v = A-kQ [2 (V^)' - ^,5{ 



Integration of these expressions over r is equivalent to integration of the left 
hand side over the invariant three dimensional volume and the final result is: 

/ dh JW^\ [A(ln V^) - STre^M] = 47rMo, 
[ dh Ji^\Aip = AnQi2-3g), 



where the notation ()44|) was used. The second integral here correspond to the 
well known Gauss law in electrodynamics. In this sense the regularization of the 
gravitational field leads to a change in the Gauss law for the electromagnetic field 
besides the above mentioned regularization. 

Another interesting result can be obtained from integration of the two differ- 
ential equations multiplied by a factor ^/gtt which gives two other invariants: 



/ dh Ji^\ [A(ln y^) - SneEM] = 47rM, 
I dh J\g(^)\A^ = 4nQ{l-2g). 
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It is obvious now that one ends up with four equations for the parameters 
M, Q Mq and g which can be solved and the latter can be expressed in terms of 
the four invariant integrals. 

9 Conclusion 

It was shown that when regular gauge conditions take place in the charged point 
source problem the solutions for both gravitational and electromagnetic fields 
are finite and have their regular behaviour at r = 0. This is a fact coming 
from the nontrivial topology of spacetime induced by the regularity of the gauge 
conditions. It is interesting to note that this topology is strongly connected to the 
definition of a physical point source described in the very beginning. It describes 
the assumption for finite structure of the point in a quite natural way. In other 
words when the condition L ^ a discussed in the beginning is no more satisfied 
the geometry has its own way of protecting the source's structure which is quite 
different from the mathematical one. This way leads to finite luminosity as well 
as finite field potentials, ip{r) and the metric components (pi^2ir), at the point 
where the source is placed. The regular behaviour of the fields is described by 
corrections to Newton's and Coulomb's laws at a classical level which take place 
only at distances defined by a new scaling parameter f. In this way inclusion 
of the gravitational field leads to a natural regularization at a classical level. 
This is a fact expected by many of the people creating GR in its beginning but 
the nontrivial geometry leading to such result is the price one has to pay. This 
geometry is essentially different from the geometry used in ^1] to describe charged 
point particles due to the difference between the Keplerian and bare mass of the 
particle. 

Another consequence of these geometrical properties is the correction in Gauss 
theorem when considering the electromagnetic field. This correction comes from 
a new parameter, recently used in description of GR problems, which is the 
gravitational mass ratio g. This is a parameter defining the geometrical behaviour 
and the new distance scale of the solution and has a clear physical meaning - it 
describes the difference between the bare and Keplerian mass of an object. It is 
a very important quantity in this problem and there are several recent attempts 
to measure its value from supernovas which will provide some evidential reasons 
for its use. 

Still after the derivation of the pure gravitational solution in |4j and the one 
for a charged point source in this article many open questions rise. The most 
important ones are the behaviour of test particles and classical fields and the cor- 
responding scattering problem, geodesic completeness and mostly the derivation 
of these solutions as a limit from finite objects. All these are subject to future 
development in this direction. 
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